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Abstract 

For the anisotropic Universe filled with massless vector field in the General Relativity frame 
we obtain bouncing solution for one of scale factors. We obtain the Universe with finite maximal 
C/} energy density, finite value of R, R^R^v, RV va P RuuaB and non-zero value of a scale factor for 

^ directions transverse to a vector field. Such a bounce can be also obtained for a massive vector 

^vq field with kinetic initial conditions, which gives isotropic low energy limit. We discuss the exis- 

tence of a bounce for a massless vector field with additional matter fields, such as cosmological 
constant or dust. We also discuss bouncing solution for massless vector field domination in n + 2 
dimensional space-time. 

£u 

^3 Introduction 

The evolution of Abelian vector fields in cosmology has been discussed in the context of generating 
isotropic inflation pQ , large scale structure [21 13 HI El [6] , dark energy and background anisotropics 
[H [7j- One shall note, that unlike scalar fields, we observe elementary vector fields in accelerator 
experiments. Thus, we expect them to appear in the early Universe, for instance as a gauge bosons. 
Since then, they are considered to be a viable alternative for scalar fields in solving problems of 
classical cosmology. As we shall show in this paper, vector fields may also result in a finiteness of 
the energy density of the Universe and generate bouncing evolution of one of scale factors in the 
GR frame. 

Our original motivation to study the evolution of a massless vector field in the GR frame was 
to find the classical limit for its evolution in Loop Quantum Cosmology [8]. However, it turned out 
that the massless vector field domination can generate a bounce already in the GR frame. Let us 
note, that bouncing solutions were found long before LQC OHDHH], for example in f(R) theories 
|12j . in the cycling Universe |13| or as a result of quantum corrections [14J. In the present paper we 
study the detailed model mentioned already in the previous paper [8J. It turns out, that the model 
has, several advantages: one does not need to introduce any modification of General Relativity 
(like in f(R) theory) or to introduce any exotic kinetic terms for fields, which would dominate the 
Universe during a bounce (like in ekpyrotic\cyclic Universe). Solutions, that we have obtained here 
may form an interesting background for particle physics models of the early Universe. 
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In this paper we will use convention 8ttG = M 2 ^ = 1. 

This paper is organised as follows: In section [T] we discuss the evolution of space-time for a 
massless vector field domination, for which we find bouncing solution of the scale factor a(t). In 
section [2] we discuss a massless vector field domination in n + 2 dimensional space-time. In section 
[3] we introduce more realistic (i.e. with isotropic low energy limit) scenarios with a bounce, such as 
the domination of a massive vector field with kinetic initial conditions or the Universe filled with a 
massless vector field, cosmological constant and/or dust. Conclusions are presented in section Hi 



1 Massless vector field domination 

1.1 Canonical vector field 

The evolution of a subdominant vector field A^ in FRW space-time has been analysed in e.g. [21 0]. 
It has been shown, that the homogeneous (background) component of a vector field lies in the 
space-like part of a four-vector. Then, one may choose one of axes to point in the direction of a 
background vector field. Let us choose it to be the z axis, which gives 

A f , = A„(t) + SA^x,t), A„ = A(t)S z fl , 5A„ = (5A t ,8A), (1.1) 

where SA^ is a vector field perturbation, which may contribute to the generation of the large scale 
structure. In this paper we will consider the Universe filled with a background vector field and a 
perfect fluid. A vector field breaks the isotropy of space, thus the metric tensor takes the form of 

gftv = Diag(l, -a{tf, -a{t) 2 , -b{t) 2 ) , (1.2) 

where a(t),b(t) are scale factors. Then the Einstein equation looks as follows 

H(H + 2U) = p = p f + p A , (1.3) 

- + l + Hn = -p± = -pf-p A , (i-4) 

a b 

2°- + H 2 = -p {{ = - Pf+PA , (1.5) 

where p/,Pf are the energy density and pressure of a perfect fluid and H = = \ are Hubble 
parameters along x, y (H) and z (H) axis. The vector field's energy density and pressure is defined 
by 

PA = ± {A 2 + m 2 A 2 ) , PA = (A 2 - rr?A 2 ) . (1.6) 

Equations of motion of a vector field and a perfect fluid looks as follow 

A + (2H - H)A + m 2 A = 0, p f + (2H + H)(p f +p f ) = . (1.7) 

First of all let us consider a massless vector field domination. This scenario has been briefly 
analysed in [8]. When the energy-stress tensor consist only of the massless vector field, then pa = 
p = p. From eq. ( 1.3|1.5 ) one obtains 

H 2 + 2HU = 2- + H 2 => - = j=> b = Eh , (1.8) 
a a b 
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where E = const > has a dimension of time. This shows, that the Universe shrinks along the z 



axis if only a < 0. Let us note, that the eq. (1.8) is also valid in the presence of additional perfect 
fluid as long as pf = —pf. The equation of motion for a massless vector field looks as follows 

p A + 4Hp A = =► p A = PAo (-Y , (1.9) 



where p Q = p(t ) and t Q is any fixed moment in the Universes history. To calculate a(t) let us 
combine eq. ( 1.5|1.9 ) which together give 

-2 



2aa + a 2 



p Q a 4 a 



(1.10) 



Let us note, that this equation is also valid in the presence of additional perfect fluid with pj = 0. 



The eq. (1.10) can be simplified to 



v6 
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(1.11) 



where pi is a constant of integration, which has a dimensionality of energy density. The x(t) is 
defined by x = a(t)/aj, where p(ai) = pj. Prom the eq. (1.11) one can see that x > 1, so pi is the 



maximal allowed energy density of the Universe [8]. So far we have made no assumptions about pi, 
but in this case we restrict ourselves to consider pj <C MS. Then GR could be the correct theory 
of gravity in the whole energy range. The exact solution of the eq. ( |1.11[ ) is 



x(t) 



2 l/3 



(2 + w 2 + y/4w 2 + w 4 ) 1 / 3 



+ 



(2 + w 2 + vW + w 4 ) 1 / 3 
2~V3 



:i.i2) 



where w 2 = |p 7 (i - a) 2 . From eq. ( |l.8|l.l2[ ) one obtains R = for any t. For t — >■ tj one also 
obtains R^ v ' R^ u — > 4p 2 and R^ ua ^ R^ va j3 — > 20pj. Thus the massless vector field domination gives 
no initial curvature singularity of the Universe Q which together with ai ^ 0, H(tj) = are features 
of bouncing solutions, like e.g. the Big Bounce. On the other hand from ( 1.8|1.11 ) one obtains 
T~L — > oo, 6—7-0 for t —> tj, which looks alike the Big Bang scenario. From now on we will set tj to 
be equal 0. When we are very close to the initial value of the energy density one obtains w <C 1, 
which gives 



x(t) ~ (1 + l -w 2 ) = (1 + l - Pl t 2 ) 



b oc t 



\ ' x I [ 1+ -pit 2 



(1.13) 



The Universe expands along the z direction for w £ (0,4). For p <ti pi one obtains w S> 1. Then 

x (t)~w 2 /*= (- A pX'\ 2 / 3 , boct-^oca- 1 / 2 , Voct. (1.14) 



The low energy limit bring us to the Kasner-like solution of diagonal Bianchi I model, with one 
direction shrinking and two expanding. Such predictions are obviously excluded by the astronomical 
observations, but in section [3] we will consider realistic extensions of this model. The evolution of 
a(w) and b(w) is shown in the fig. [Tj 

J The energy conditions for the vector field domination require p > 0, p > and p > \p\. The massless vector field 
satisfies them, since p = p = A 2 /2b 2 > 0. 
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Figure 1: Left and right panels show the evolution of a(w)/aj and b(w) / '(aiE) respectively for 
massless vector field domination. Let us note, that 6(0) = 0. In both panels solid blue lines and 
dashed red lines represent the evolution of scale factors and its low energy limit respectively. 



1.2 Vector field non-minimally coupled to gravity 



Let us investigate the case of a vector field with non-minimal coupling to gravity. This kind of 
coupling has been introduced to provide the slow-roll evolution during inflation for the physical 
vector field defined by V(t) = A(t)/b [HI21CE5]- In such a case a vector field may survive inflation 
and dominate the Universe after reheating. The action for a non-minimally coupled background 
vector field looks as follows 1161 



S A 



(Fx 



1 



(m 2 + £R)£ 



S = Sa + Sj, 



11.15) 



where g = det(g^ u )- Let us assume, that the perfect fluid is minimally coupled to gravity. Then one 
obtains the effective Einstein equation in the form 



H(H + 2H) = p f + ^ + ^(H + 2/7 ), 1,1 + (nr + 2< 



a b 



2- + H z 
a 



-P± 



2b 2 



(1-40 



A 2 



-2U 
2£ 



A 2 



2b 2 b 2 



2UH + i/ 2 )) 2&2 , 
A((3H- H)A -A] + 



m 2 + 2( [m 2 



UH + H 2 -U + H 



-p/ + (l + 4£) 



A 2 2^ 
2b^~¥ 



A [(AU- 2H)A - A 



A 2 ^ 
2b 2 



m 2 + 2£ 



-2V 2 + 87iH + 3H 2 + m + 2H 



J1.16) 



:i.l7) 



11.18) 



The bouncing solution appears only if initially energy of a vector field's kinetic term dominates, 
since we w ant to avoid bit) dependence of the parallel pressure and the energy density. From eq. 
( 1.16[l.l8 ) for A 2 3> max{m 2 A 2 , RA 2 } one obtains approximate equations of motions of scale factors 



x 



X 
X 



b x 
b x 



2ipix~ 



J1.19) 
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Thus, from the eq. ( 1.11 ) one recovers (1.12 ) as a solution for x(w) around bounce, but the definition 



of w variable is changed to w = |p/(l + 4£)i . We have obtained bouncing solution for x with 
slightly changed effective time variable, but we need to assume, that 1 + 4^ > to keep w real. Let 



us note, that the non-minimal coupling violates the eq. (1.8). To find the evolution of b for energies 



close to pi let us consider b = b(x). Then from the eq. (1.19) one obtains 



1 _ l+2€ 

b oc (x - l)2(i+40a: ->• for x -> 1 , (1.20) 
so one recovers all features of bounce generated by a minimally coupled massless vector field. 

1.3 Vector field with non-canonical kinetic term 

Another way to obtain the slow-roll evolution of pa during inflation is to consider a vector field with 
a non-canonical kinetic term. Such a field may dominate the Universe after inflaton and become 
the vector curvaton (5j E] . Let us define its action as 



<l l .'\/-fi ( \R - f\F^ + im 2 .*".!, ) . (1.21) 



where / = f(a, b),m = m(a, b). Time dependence of the mass term is necessary to obtain a constant 
value of M = m/y/f, which is a mass term of a canonically normalised vector field U = xffA/b. 
Einstein equations, together with the equation of motion of U, look as follows 

H(H + 2H) = p f + p A = Pf + ^(fA 2 + m 2 A 2 ) , (1.22) 
lA + nH --^-^ + ^(f + H)-f^ + ^) ■ (1-23) 

2 l + H2 = -n = -pf~w^ f - bf > b) + w^ m2 ~ b ^ ' (L24) 

U +{2H + H)U + (m 2 + 2UH + U + ^-^_^2H-H)j) U = 0. (1.25) 



To obtain cancellation of all time dependent mass terms in the eq. (1.25) one needs / oc b , m oc b. 
For / oc a -4 , m oc a~ 2 the time dependent part of the mass term is equal to 2H+7-1, so it is negligible 
during inflation. Let us note, that for both sets of /, m functions one obtains p± = pu , so T% is fully 
isotropic [^J The continuity equation for a vector field is of the form of 



PA + {2H+H)(p A +p A ) = for / oc b J , PA + {2H+H) {p A -p A ) = for / oc a" 4 . (1.26) 

One can see, that / oc b 2 gives T% identical with a scalar field domination scenario. 

On the other hand, for / oc a -4 one obtains surprising results. For the kinetic term domination 
one obtains p A — PA, which means that pA *C Hp a- This case is similar to the slow-roll phase of a 
scalar field. The potential term domination brings us to pa — —PA, so pa — —6HpA, which is the 



2 In \5\ [B] Dimopoulos et al consider m oc a and / oc a 4 or / oc a 2 to obtain flat power spectra of initial 
inhomogeneities . 
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evolution similar to the massless scalar field scenario. In general, one can see, that the behaviour of 
kinetic and potential terms has been swapped, comparing to a massive scalar field scenario. 

For example, for a massless vector field (or massive vector field with kinetic initial conditions) 
with / oc o -4 one obtains pa = => pa = const. Then, from the eq. (1.24), one obtains 



2- + H' 
a 



PA 



a 
ai 




(1.27) 



We have obtained a bounce for a = aj. From now on we shall replace pa by pj, since pa does not 



need to be constant in the massive case. The solution of the eq. (1.27) looks as follows 



x(t) = cosh 2/3 



3pj (t - tj) 



cos 



h 2/3 



W 

7s. 



(1.28) 



where w = 3^/~pJ(t — tj)/2. For a massless vector field, from eq. qi.22|1.24| ) one obtains b oc a oc 
x, so we have found analytical solutions for both scale factors. For such an evolution of scale 
factors one obtains no initial curvature singularity, since R — >■ — 51/7//4, R^ V R^ U — > 1413/)|/32, and 
R^ va P R^vap — > 549/7^/16 for x — > 1. Let us note, that for big values of w one obtains 



a(t) oc b(t) oc eVw/ 5 * . 



(1.29) 



This exponential expansion is identical with the one of the FRW Universe filled with a cosmological 
constant. Thus, such a vector field can be a good candidate to play the role of Dark Energy. Let 
us note, that the Universe becomes isotropic very fast, even without a vector field's mass term|^] 



Let us note, that for a = b the a(t) from the eq. (1.28) does not satisfy 3H = pi = const. Thus 



anisotropic initial conditions for scale factors are necessary to obtain bouncing solution for a(t). 

On the other hand, one can ask about the naturalness of a kinetic term's domination at early 
times. For a scalar field it is natural from the point of view of the phase space to consider the kinetic 
term domination around a bounce. Thus, at early times, due to analogy with scalar fields, one shall 
expect the domination of a vector field's potential term for / oc a -4 . This means, that the bouncing 
solution is strongly fine tuned if M. ^ 0. Let us also note, that the bouncing solution described by 



the eq. (1.28) may be also obtained with a cosmological constant. In particular, it may be obtained 



for a scalar field with potential initial conditions. Thus, such a bouncing solution is not a unique 
feature of vector fields. 

The evolution of scale factors and Hubble parameters for a massless vector field with non- 
canonical kinetic term is shown in the fig. [2} 



2 Vector field in n+2 dimentions 

As we have proven, a massless vector field generates a bounce for a(t) and gives the Kasner-like 
evolution at the low energy limit. Let us consider more generic case of a massless vector field 
domination in n + 2 dimensions, where n is a number of dimensions with a rotation symmetry 
Let us assume, that the vector field points x n+ \ axis, which gives the metric tensor of the form of 

3 Or with the negligible influence of the mass term on the evolution of the Universe. 
4 This means, that the 4 dimensional space-time is described by n = 2. 
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Figure 2: The left panel shows the evolution of scale factors for the domination of a massless vector 
field with non canonical kinetic term and / oc a -4 . The right panel presents the evolution of Hubble 
parameters in the same model. One can see, that after w ~ tJ^H ~ 5 the Universe becomes isotropic 
and enters the era of exponential expansion. 



gik = Diag (l, —a 2 , —a 2 , . . . , —a 2 , — ft 2 ), where l,k = 0, 1, . . . ,n + 1. Then from Einstein equations 
one obtains 

1 "1 CL I) 

-n(n - 1)H 2 + nHU = n- + -n(n - l)H 2 ^- = -^b = Ea. (2.1) 

2 a 2 do 



It is worth to note, that the eq. (2.1) is valid for any n. From the continuity equation one obtains 
p oc a~ 2n , so from G n ^+i = Pa one obtains the equation of motion for a scale factor a(t) 



1. . 2 

xx + -(n — l)x 



Pl x 2(n-1) 

n 



(2.2) 



where x{t) = a(t)/aj, ai = a(tj) and ti is the moment of a bounce. Let us assume, that n / 1, 
Then, from the eq. (2.2) one finds 



x = x 



l-n 



2 PI 



n{n — 1) 



\Zx n ~ x - 1 



(2.3) 



The eq. (2.3) has the solution of the form of 



%/2(w- l)ry 1+n)/ 



' 2-^1 



1 + n 1 n — 3 



„i-n 



2(1 -n) } 2' 2(n- 1)'' 



(2.4) 



where 2-^1 is a hypergeometric function. Analogously to the 4 dimensional case a bounce appears 
for i = tj. The is a constant of integration and the maximal energy density of the Universe. The 
eq. (2.4) can be simplified in two limits. Let us consider the evolution of scale factors around a 
bounce (which means, that x — 1 <C 1) and in the low energy limit (for which j;> 1). Then from 
the eq. (|2.3l) one obtains 



x~l + — (t-t/) 2 for 
2n 



x ~ for x S> 1 



(2.5) 
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Those simplified solutions are valid for any n. The n = 1 case gives the equation of motion of the 
form of 

x = ^2 Pl \n[x(t-t I )\ . (2.6) 
Thus, for x = 1 one obtains x = 0. In such a case there is no analytical solution for x(t). 



The initial curvature singularity of the Universe 

To discuss the issue of initial singularity in n+2 dimensions let us express R, R^R^ and R^ ua P : R^ V0l j3 
us a function of x. From eq. ( 2.1|2.3 ) one obtains 

R = l( n -2)x- 2n Pl , (2.7) 
n 

R^R = 7 4, (2n 2 - 3n + 2) x~ An pj , (2.8) 

R^RuvaB = ^ ^ 9 ( n 2 (n + l)x 2(n - 1} - 4n 2 (2n - l)*"" 1 + 2 (8n 3 - 12n 2 + 7n - 2)) .(2.9) 
(n — l)rr V / 

In the limit of a bounce, i.e. for x — > 1, one obtains 

R -»■ 2 (n-2)p/ , fl""^ ^ 4 ( 2n2 " 3n + 2) , R^R^ afi -)■ 4 ( 9n2 ~ 10n + 4 ) P? ■ 

(2.10) 

Let us note, that for n ^ all invariants mentioned above are finite around a bounce. Thus, for any 
n 7^ one obtains no initial curvature singularity of the Universe. 



3 Isotropic models with classical bounce 
3.1 Massive vector field domination 

The fact, that a massive vector field with kinetic initial conditions can give bouncing solution for 
a(t) has been proven in [8j. During a bounce the kinetic term dominates over the potential one. 
Nevertheless, A(t) grows with time, which leads to vector field's oscillations. During that period the 
background anisotropy {H — H)/(H + T-L) decreases like t~ l and the Universe becomes isotropic. To 
find more about the evolution of space-time in such a model see [Sj • 



3.2 Massless vector field and cosmological constant 



Analytical solutions of the Einstein equation for the massless vector field domination scenario may 
be obtained also for a massless vector field with a perfect fluid as long as pf is not a function of b(t). 
Then both sides of the eq. (1.5) are functions of a(t) and its time derivatives. The simplest case is 

-A 



a cosmological constant, which gives pf 
tj, one can simplify the eq. (1.5) to 



-Pf = PA 



const. Then, with a certain choose of 



x(x) 



1 



t) x + y x PAI - 



PAI 



(3.1) 
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From the eq. (3.1 ) one obtains x — > for x — > 1, so aj is a minimal value of a(t). Thus, one obtains a 
bounce similar to the one from a massless vector field domination. The pi = p\ + pai is the maximal 
value of the energy density, since pa = const and pa = pai%~ 4 - Let us note, that a cosmological 
constant, together with a massless vector field, has along the z axis the effective equation of state 
p = —p. Then, form eq. ( 1.3|1.5 ) one obtains b = Ed, so the eq. (1.8) is still valid. Thus, there is 
a clear analogy between cases with and without a cosmological constant. Simply, one expects the 
cosmological constant to dominate the Universe for energies much lower than pj, so it shall not have 
significant influence on the evolution of space-time around the bounce. To find the low energy limit 
of that model let us assume, that i> 1 and pa > 0. Then one finds 



x(t) oc sinh 2 / 3 ( ) -> e^~ 3t 



b oc x 



(3.2) 



so for x 3> 1 one obtains isotropic Universe even without a vector field's oscillations. 



For the Universe filled with massive vector field with kinetic initial conditions and with a cos- 
mological constant one obtains the following low energy limit for the evolution of space-time: When 
a vector field starts to oscillate the Universe becomes isotropic. Vector field's oscillations generate 
particles and radiation, so in the presence of A one obtains the observed Universe. This realistic low 
energy limit brings us a chance to consider a situation, in which there is no need for any theory of 
gravity wider than GR. Let us also note, that a cosmological constant may come from the massless 



vector field with non-canonical kinetic term (see section 1.3). Thus, one can obtain a bounce with 
a good low energy evolution of the Universe using only vector fields. 



3.3 Massless vector field and dust 



The eq. (1.12), which is the analytical solution for a(t) obtained for the massless vector field 



domination, is also valid for the Universe filled with a massless vector field and dust, since dust 



does not contributes to py. Thus, the eq. (1.11) remains unchanged. On the other hand, from the 
continuity equation, one obtains pa ~ a -4 , p m ~ a _2 6 _1 , so the energy density is a function of both 
scale factors. This means, that even if one defines pai such as V a pa < Pai, there may be no p m i 
such that V a .(, p m < Pmi, so the total energy density does not need to be finite at all times. Let us 
note, that the eq. (1.8) is not valid, since from eq. ( 1.3|1,5 ) one obtains 

HH - - = \ Pm ± . (3.3) 
a 2 

This equation may be written as a equation of motion of b(t). Let us define x = a/aj, y = b/b Q , 
%o = a Q /ai where aj = a(ti), b Q = b(t ), t Q ^ ti is a certain moment and tj is the moment of a 
bounce of a(t). The t Q needs to be different than ti since we do not know, if p m is finite at t = tj. 



From the eq. (3.3) one obtains 



xy -xy = p rno 



2x' 



From eq. fll.ll|3.4 ) one finds the evolution of y as a function of 



x 



y(x) 



Pfo 



3pAi x 



-(x 2 + 4x 



+ C-v / aT 
x 



(3.4) 



(3.5) 
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where C is the constant of integration. Thus, for x — > 1 one obtains (x 2 + 4x — 8)/x — > —3 and 
yx — 1/x —7- 0, which means, that for x close enough to 1 one obtains y < 0. This means, that the 
minimal allowed value of x need to be larger that 1. Let us define x p such, as y(x p ) = 0. Let us note, 
that H(x p ) > 0, since x p > 1, which together with H(x p ) = oo gives H(H + 2H){x = x p ) = oo. This 
seems to be natural, since p m oc (x 2 y)^ 1 =^ p m (x p ) = oo. Thus, by adding dust to a massless vector 
field one obtains bounce of a(t), which happens in the forbidden region of y < 0. The beginning 
of the Universe appears for y = 0, for which the energy density is infinite. One cannot extend 
the evolution of the Universe for t < t p , so the considered model has all features of the Big Bang. 
The only difference is, that the beginning of space-time is not a point, but a circle described by 
y = 0,x = x p . 

Similar results can be found for the Universe filled with massless vector field, cosmological 
constant and dust. There for x — > 1 one obtains y — > —p mo x 2 /(pAi + A), so for A < and pa < —A 
one obtains y > for x < x p . Thus, one can solve the problem of negative y at the moment of 
bounce. However, this generates another problem, since y becomes negative for x > x p . In such a 
case the Universe would start with a bounce to immediately end up in the Big Crunch. This means, 
that y{x) may be positive for x < x p or for x > x p , but one cannot obtain y > at all times. 

4 Conclusions 

In the initial section of this paper we have shown, that a massless vector field domination generates 
bouncing evolution of space-time in a plane transverse to a vector field's direction. The bouncing 
solution us has many features of other bounces known in cosmology, such as a non-zero initial value 
of the scale factor a(t), lack of initial curvature singularity and finite initial (maximal) energy den- 
sity. Along the vector field's direction one obtains the evolution of the scale factor and the Hubble 
parameter similar to the Big Bang model, i.e. — ?- and T~L — > — oo for t — > tj. The Kasner-like 
solution with a ~ t 2 / 3 , b ~ t- 1 / 3 is the low energy limit of the evolution of space-time. 

In section [2] we have considered bouncing solution for a massless vector field domination in the 
general n + 2 dimensional case, where n is a number of dimensions with isotropic evolution. We 
have found the analytical solution for the scale factor a(i), as well as its simple form in the high and 
low energy limit. We have also proven, that R, R^R^ and R^ VOL ^ R^ap do not diverge for t — > tj 
for all 

In section [3j to obtain isotropic low energy limit of the evolution of space-time, we have discussed 
several extensions of the massless vector field domination scenario. All of them obtain bounce at 
(x,y) plane for t = tj. First of all we have considered massive vector field with kinetic initial 
conditions. I such a case the Universe has bouncing solution for p ~ pj and isotropic low energy 
limit, which is generated by the vector field's oscillations. Secondlyl we have considered the Universe 
filled with a massless vector field and a cosmological constant. In this case a cosmological constant 
contributes to pi and generates isotropic expansion for a 3> aj. The massive vector field with 
cosmological constant may give realistic low energy limit with particles and radiation originating 
from the and with a Dark Energy coming from A. We have also considered the Universe filled 
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with massless vector field and dust, for which the bounce appears in the forbidden region b < 0. The 
true beginning of space-time appears at b = 0, H = oo, a > ai,H 0, so the initial energy density 
diverges and one obtains all features of the Big Bang. Similar results are given by a massless vector 
field, dust and cosmological constant domination. 

Let us also note, that the presented model has several advantages comparing to bounces from 
e.g. f(R) theories or ekpyrotic\cyclic Universe. GR remains the only theory of gravity needed and 
the matter sector has the standard form known from particle physics. 
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